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is found to exhibit a second-order transition with logarithmic singularities at criticality. 
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I. INTRODUCTION 

The Ising model in a nonzero magnetic field is a 
well-known unsolved problem in statistical mechanics. 
In 1960, Fisher pj produced a remarkable solution of 
a super-exchange antifcrromagnetic Ising model in the 
presence of a nonzero field. The Fisher model is defined 
on a decorated square lattice where there is an external 
magnetic field applied to decorating spins which interact 
via super-exchange interactions. 

In this paper we consider similar super-exchange mod- 
els on the kagome lattice. The kagome lattice has been of 
interest in recent years in the role it plays in the high-T c 
superconductivity. It has been known that special cases 
of the Kagome Ising model are soluble in the presence 
of a magnetic field 0, Eg ■ The structure of the ferrimag- 
net SrCrgGa40ig, for example, is found to consist of 2D 
spinel (kagome) slabs with magnetic spins residing at 
1/3 of the lattice sites. Thus it is of interest to consider 
models with similar structures. Azaria and Giacomini 
have extended the Fisher model by considering the 
kagome lattice shown in Fig. [J a), which reduces to the 
Fisher model upon setting K\ = 0. They obtained its 
partition function from which the phase diagram is de- 
duced. But there has been no detailed discussion of its 
thermodynamic properties. Here, we consider yet an- 
other extension of the Fisher model to a kagome lattice 
as shown in Fig. ^b). We show that the partition func- 
tion of this model is identical to the Azaria and Giaco- 
mini solution p| , although there appears to be no direct 
mapping between the two models. In addition, we also 
present detailed analyses of the thermodynamics of the 
solution. 

Consider an Ising model shown in Fig. ^b) with an 
interaction energy 

— J\d\U% + J((72<Jz — <73<7i) 

around every triangle formed by spins <7i, 02,03. Intro- 
duce reduced interactions K — ft 'J, K\ — j3J\ , where 
[3 = 1/kT, such that K,K\ > indicate ferromag- 
netic interactions. In addition, there is an external 
magnetic field H applied to 2/3 of the lattice sites de- 
noted by solid circles. We denote the reduced field by 



L = H/kT. As a result, the magnetic spins interact 
with an super-exchange interaction via intermediate non- 
magnetic spins. It is clear that there is no loss of gener- 
ality to restrict considerations to 

H,J>0, or L,K>0. 

For K\ = the model reduces to the Fisher model 
for the square lattice. For K\ 7^ the present model 
is more general and differs from the Fisher model in a 
fundamental way as described in Sec. 2 below. 

Our main result is a closed-form solution of the par- 
tition function and detailed analyses of thermodynamic 
and magnetic properties. In Sec. 2 we deduce the so- 
lution using a combination of star-triangle and decima- 
tion transformations. The phase diagram is analyzed in 
Sec. 3. In Sees. 4-6 the internal energy, specific heat, 
magnetization, and susceptibility are analyzed using a 
free-fermion model formulation [6j. 



II. THE PARTITION FUNCTION 

Denote the partition function of the kagome Ising lat- 
tice in Fig. 1(b) by Zkg(K, Ki, L). Our main result is 
an equivalence of Zyx^{K,K\,L) to the partition func- 
tion of an honeycomb Ising model in zero field, a result 
which renders the model soluble. 

This equivalence is established by effecting a sequence 
of spin transformations. First, we carry out a star- 
triangle transformation for every triangular face as shown 
in the first line in Fig. 2. This converts the kagome lattice 
to a decorated honeycomb lattice. Next the decorating 
spins are decimated as shown in the second and third 
lines in Fig. 2, and the lattice is reduced to that of an 
honeycomb. The crux of matter is that the external fields 
L\ and —L\ induced in the second step cancel out at the 
end. As a result, the final honeycomb lattice has no ex- 
ternal field and is soluble. It has reduced interactions 
R, R and R\ in the three principal directions. 

The transformations shown in Fig. 2 are standard 0, 

. For the transformation in the first line in Fig. 2, we 



2 




(a) (b) (c) 

FIG. 1: (a) and (b) are two soluble Kagome lattice Ising models in a magnetic field. Solid circles denote magnetic spins and 
open circles non-magnetic spins. Both models are mapped to an honeycomb Ising lattice without a field shown in (c). F N is 
an overall factor. 
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FIG. 2: Transformations used in deducing the mapping shown 
in Fig. 1. Fi, i = 1, 2, 3 are constants. 



have 



-2K-K t 



2K — K\ 



2F X cosh(Ti + 2r) 
2F l cosh(ri - 2r) 



e Kl = 2FiCOshri 



(1) 

(2) 
(3) 



from which we can solve T and T± in terms of K and K\. 
Divide the sum and difference of and J2Jl by 10, we 
obtain, respectively, 



tanhTi 
We also obtain 



cosh2r = e-^cosh^if 
e 



~ 2Kl sinh 2K 



sinh 2r 



(4) 
(5) 



F* = [16 cosh 2 TiCOshiTi + 2r)cosh(ri - 2r)] \ (6) 

Now K > and Ki real, so (QJ and JSJ show that there 
are two regimes 

e~ 2K 1 cosh 2K > 1 regime I 

e~ 2K 1 cosh 2K < 1 regime II (7) 

and r and 1\ are real in regime I and pure imaginary in 
regime II. 



For the transformation in the second line in Fig. 2, we 
have 



2cosh2ri = F 2 e 

V 2 



Hi 



2 = F,e" fll 



from which we obtain 



F-2 = 4cosh2ri 



(8) 



(9) 



e 2fll = cosh2r 1 = £ _ 4Ki _ 1 (10) 

where the last step is obtained after making use of J3J. 
Equation (|10|l shows that i?i is real, and R\ > in regime 
I and i?i < in regime II. 

For the transformation in the third line in Fig. 2, we 
have 



2cosh( J L + 2r) = F 3 e- R+2Ll 
2cosh(L-2r) = F 3 e- R - 2Ll 
2 cosh L = F 3 e R 



(11) 
(12) 

(13) 



from which we obtain 



_ 4R _ cosh(£ + 2r) cosh(L - 2r) _ 
cosh L 
e- 4Kl cosh 2 2K + sinh 2 L 

cosh 2 L 
4Ll = cosh(£ + 2r) 
cosh(i-2r) 

Fi = 16cosh 2 Lcosh(£ + 2r)cosh(L-2r). (16) 



sinh 2 2r 
cosh 2 L 

(14) 
(15) 



Thus, R is always real, and we have 

R < 0,Ri > 0, T, ri, Li = real, regime I 

R > 0,Ri < 0, T, Ti,Li = imaginery, regime 11(17) 

This says that R and R± have opposite signs. In the 
solution obtained by Azaria and Giacomini |5j], the final 
equivalent honeycomb Ising model is identical with ours 
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except the replacement of R by —R. However, as we 
shall seen l|23[) below, the negations of R and/or R x do 
not affect the solution. Therefore, the two solutions are 
identical. 

Combining ©, ifR)|) . and after some reduction, we 
obtain 



F = 



4 cosh 2 L(e 4Kl sinh 2 L + cosh 2 2K) 
cosh 2 i?i 



(18) 



Now the field L\ at the decorating sites is canceled and 
we arrive at the final equivalence 



Z KG (K,K 1 ,L)=F N Z HC [R(L),Ri 



(19) 



where Znc[R(L), Ri] is the partition function of the hon- 
eycomb Ising lattice. Here, N is the number of triangular 
faces (or the number of magnetic spins) of the kagome lat- 
tice, which is also the number of honeycomb lattice sites. 
Note that the L dependence is only in R. 

The transformation 11911 reduces to that of Fisher's for 
the square lattice after setting K\ — [R\ = oo). But 
our model differs from the Fisher model in a crucial as- 
pect: The Fisher model results in a square lattice whose 
sites are the decorating sites of the original lattice, in 
the present model all sites of the original kagome lattice 
disappear at the end. Thus the ordering of the honey- 
comb lattice bears no direct relationship to that of the 
kagome lattice. In addition, while there is no frustrated 
plaquettes in the Fisher model, in the K\ > model all 
triangular faces are frustrated. 

Introduce the per magnetic spin free energy 



/ KG (^i,£) 



= lim N-HnZKciK^^L) 

iV— >oc 



and the per-site honeycomb lattice free energy 

fuc(R,Ri) = lim N-HnZnciR.Ri)- (20) 



jV- 



Taking the N — > oo limit of l|19f) and making use of the 
explicit expression of /hc given in H, EJ , we obtain for 
the kagome Ising model, 

f KG (K,K u L) = In F + fuc{R,Ri) 



= lnF+ -ln2 
4 



1 



167T 2 Jo 



2;r 



lnS(M) (21) 



where 



= cosh 2i?i cosh 2 2R+1- sinh 2 2R cos( 
— sinh 2Ri sinh 2i?(cos 9 + cos </>) . 



(22) 



Note that the negation of either R or Ri corresponds to 
changing 6 ^ tt — 9, <f> — * tt — 4> in 11221) which does not 
change the free energy /kg- 



As a check, we recover the Fisher solution 0] upon 
setting K\ = or R\ — > oo: 

/ KG (K, 0, L) = | In 2 + i ln[cosh 2 L(sinh 2 L + cosh 2 2K)] 



1 



16tt 2 



In 



cosh 2 2R 
- sinh 2i?(cos 9 + cos (f>) (23) 



D -AR _ 



(sinh 2 L + cosh 2 2K) / cosh 2 L. 



with 

The internal energy and magnetization per magnetic 
spin are, respectively, 



U(K, K U L) 



d_ 

df3 



f KG (K,K u L), 



(24) 



M(K,K U L) - —fgG^K^L). (25) 

The specific heat and susceptibility are further computed 

as 

C H (K,K U L) = —UtfCKuL), (26) 
x(K,Kx,L) = -^M(K,K U L). (27) 



III. THE PHASE DIAGRAM 

The integral in the free energy /kg is precisely of the 
form of that of the free-fermion model discussed by Fan 
and Wu for which the critical condition is (for RRi < 
0) 



S(7T,7T) = 0. 



(28) 



Alternately, the critical point for an anisotropic honey- 
comb Ising lattice has been given in 9j as C1C2C3 + 
1 = SiS 2 + S2S3 + #3 Si where Q = cosh2R l ,S l = 
s'mh\2Ri\,i = 1,2,3. It is also given by the expression 
uji = ui2 + + LU4 where the oj's are given in the next 
section. Using any of these expressions, we obtain the 
critical condition 



smh[2R(L)} = -coth-Ri. 
Explicitly, lt2TJ|l reads 



cosh L 



(29) 



1 + tanli Ri - 1 
e- 4Kl cosh 2 2K - l) , T = T C {H$Q) 



which can be realized only in regime I and K\ < K , or 
7 < 1 (see below). Note that for K 1 = lf3*0"|l reduces to 
the Fisher expression 



coshL = J(V2 - l)/2 smh2it" 



(31) 
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It can be verified that we have 

|sinh2iJ(i)| < | cothiii|, T>T C (H). (32) 

Introduce parameters 



a = L/2K = H/2J > 
7 = J 1 /J = K 1 /K 



(33) 



such that 7 > indicates J\ is ferromagnetic, and con- 
sider the phase diagram (|30l) in the {a, 1/K} plane, 
where 1/K is the temperature. The phase diagram is 
plotted in Fig. 3 for different fixed values of 7. 




1.5 2 
Temperature 1/K 

FIG. 3: Phase boundary in the a-T plane for fixed values of 
7- 

At low temperatures the phase boundary behaves as 

(34) 



1-7-^^2(^2 + 1) 



with an initial slope independent of 7. Solving (|30|) for 7 
at H = 0, we obtain 



7 



1 

2K r 



2 ln(sinh 2K C ) - ln(2 cosh 2K C ) 



(35) 



We see that T c (0) decreases with 7, reaching T c (0) = 
at 7 = 1. It follows that there is no transition when 



7 > 1, or K\ > K 



(36) 



as aforementioned. This can be understood physically 
since, in the limit of 7 — > 00, the non-magnetic spins 
can assume values ±1 randomly and there is no ordered 
state. 

For small H l|30|) is expanded as 



T C {H) =T c (0)-CH 2 , 
where C is a constant given by 

32(16 - 24< 4 + < 8 ) 16(-4 7 + 4i + t 5 ) 



(37) 



C'- 1 = kJ 



t 9 (4 + 4t 2 -t 4 ) 



16 



with t = tanh2J/A:T c (0). 

To determine the nature of regimes in the phase dia- 
gram, we consider the ground state. At zero temperature 
the kagome lattice can assume two ordered states: 

(i) All magnetic spins aligned in the same direction 
with the total energy 



-N{2a + ~f)J 



(38) 



independent of the spin directions of the non-magnetic 
spins. Since the N/2 non-magnetic spins can point in any 
direction, the entropy of this ground state is |iVfcln2. 

(ii) The magnetic spins connected by interactions — J 
assume the value +1 (or —1) and those connected by 
interactions +J assume the value —1 (or +1), while all 
non-magnetic spins are fixed at +1 (or —1). This forms a 
super-exchange antiferromagnetic state with the energy 



U 2 = -N{2 - 7) J, 



(39) 



and the ground state is two-fold degenerate with zero 
entropy in the thermodynamic limit. 

Comparing U\ with Ui-, for a > 1 — 7 which is the 
regime to the "exterior" of the critical curve (|30() . we 
have U\ < U2 indicating the system is paramagnetic. 
For a < 1 — 7 which is the regime enclosed by the critical 
curve, we have Ui > U2 and the system assumes a super- 
exchange antiferromagnetic ordering. 



IV. THERMODYNAMIC PROPERTIES 

To analyze thermodynamic and magnetic properties 
given by (|2"lf> - (|2"7Jl . we make use of results of [g. Writing 
(|2~2"|l in the form of Eq. (16) of 0, we have 

E(6,<f>) = 2a + 2b cos 9 + 2c cos </>+ 2d cos(0 + cj>), (40) 

with 

1 



a 



cosh2i?i cosh^ 2R + 1 = -(wf + w| + u l + 



b = — sinh 2i?i sinh 2R = LO1UJ3 — 

c = — sinh 2i?i sinh 2i? = — CJ12W3, 

d = — sinh 2 2i? = — W1W2, 



ui = cosh(2R-R l ) 
cu 2 = cosh(2i? + i?i) 
0J3 = uja = cosh R\ . 



(41) 



(42) 



Using l|42JI we compute the parameters {a;, y, z} intro- 
duced in [6j as 



X = LU1LJ4 — Ll>2^3 



sinh 2R sinh 2Ri 



y = i — lo\ — uj\ + Lof) = — - sinh AR sinh 2Ri 



LU1OJ4 + (^2^3 = 2 cosh 2R cosh R\. 



(43) 



5 



Here we have reversed the sign of x from [6( to make x > 
0. This is permitted since only x 2 appears in the ensuing 
discussions. The critical condition l|29|l is equivalent to 
y = z and we have the regimes 

y> z> x>0, T < T C (H) 
z>y>x>0, T> T C (H). (44) 

It was established in |(| that derivatives of the free 
energy are best computed by first carrying out one-fold 
integration. Adopting notations in |6j, we obtain after 
differentiating H21|l the expression 

[f KG (K, K 1 , L)}' = C + C x h + C 2 I 2 + C 3 I 3 , (45) 

where the prime denotes the derivative with respective 
to some variable such as T and H, 

Co = 



where 



and 



Ci = 

c 2 - 
c 3 ■- 

h 
h 
h 



F' 


6' 


~F 


+ W 


F' 


d! 


~F 


+ 4d< 


f 

a - 


a fb> 


~2\b 




b' d' 


K 


~b~~d 



b 2 -d 2 (b' d' 



Ad 



1/6' d'\b 2 
2\ ~b ~ ' 



2d 



2nd 
b 2 + d 2 



- 1 



(46) 



1 

8^ 
1 

8^ 



d<t>[Q(<t>)] 



-1/2 



d<f> cos (j)[Q{(f))] 



-1/2 



1 [ 2V 

— / d<f>(l + cj cos <ft) 1 

87T Jo 



-1/2 



(47) 



with 



uj = -2bd/(b 2 + d 2 ) > 
Q{<f>) = x 2 (coscj) — yz/x 2 ) 2 

+ (x 2 -y 2 )(z 2 -x 2 )/x 2 . 



(48) 



Carrying out the integrations in l|47[). we obtain explic- 
itly for T>T C {H) (z > y > x > 0) the result 



h = 
h = 
I: 



2w(z 2 -x 2 ) 1 / 
1 



2Tryz(z 2 - x 2 ) 1 / 2 
1 



z 2 K{k) + (y 2 - z 2 )n(r, k) 
yK(k) 

'n(s,fc) (49) 



2w(y + Luz)(z 2 -x 2 y/ 2 

| u{z 2 -y 2 ) 
z + wy 



k 2 
r 
s 



(y 2 -x 2 )/(z 2 -x 2 ), 
{uz + yf/iz + ujy) 2 



(50) 



and K and II are elliptical integrals of the first and third 
kinds, respectively, 



K{k) 
IT(r, k) 



*' 2 da 



1 1 \/\ — k 2 sm~ a 



tt/2 



1-2 ci^ 2 ( 

da 



o (1 — r sin 2 a)yl — k 2 sin 2 a 



(51) 



For T < T C (H) (y > z > x > 0) we obtain the same 
result with y and z interchanged. 

The apparent non-analyticity of Co at \Ri\ — \R\ as 
indicated in l|46() is spurious. It can be shown that the 
combination of Co + C3/3 is always analytic. 

At the transition temperature y — z (or k = 1, r = 1) 
both K and II diverge as In \T - T C (H)\. Applying gSJ 
to the internal energy (|24|l where the derivatives are with 
respect to T, we obtain 

U(H) = U C (H) + const |T - T C (H)\ In \T - T C (H)\, 

T^T C (H). (52) 

A further derivative of U(H) as given by l|26l) gives the 
specific heat Cr- Thus, the energy is continuous at 
T C (H) while the specific heat diverges logarithmically. 
These findings, which are the same as those found in the 
Fisher model, indicate the occurrence of a second order 
transition along the phase boundary in Fig. 3. We plot 
U(H) and Ch respectively in Figs. 4 and 5. For com- 
pleteness, we derive the explicit expression for U C (H) in 
the next section. 

In the case of zero magnetic field H = a = and 
7 = — 1, the internal energy assumes a simple expression 
given by 



U(0)/J= -1 + (1 + tanh2A')|l - ^tanhi? 

- 1 + 2(cosh 3 2i? - 3cosh2i?- 2)/i coth2i? 

with er 2R = e 2Rl = {e AK + l)/2. 

V. DERIVATION OF U C (H) 
Define an auxiliary variable 



j(53) 



<Pa(k) = sin 



1 — a 
1-k 2 



> 



6 




0.5 1 1.5 2 2.5 

Temperature 1/K 

FIG. 4: Energy U versus temperature T at fixed magnetic 
field a for 7 = — 1. The broken line indicates the critical 
energy U c . 




FIG. 5: Specific heat Ch as a function of temperature T at 
fixed magnetic field a for 7 = — 1. 



which gives, for fc 2 , r, s given in l|5()jl. 



(f r (k) = sin y 1 — 



• -1 (1 ~ C0 2 )(z 2 - X 2 ) 

ip s (k) = sin W 



(z + wy) 2 



V3 S (1) = sin 1 



1 - w 

l+UJ 



, T = T c (ff)(54) 



where the last line holds at the critical temperature fc = 1 
or y = z. 

For k 2 < r < 1 we have [111 



II(r,fc) = K{k) + 



2 V (l-r)(r-fc 2 ) 



1 - A (ip ri k) 



where Ao is Heuman's Lambda function 
A (tp, k) = - \K(k)EUp, fc') - \K[k) - E(k)]F(ip, fc') 

7T L 

(55) 

with fc' = yl — fc 2 and 



(56) 



£(<M) = / (1 - fc 2 sin 2 6») 1/2 d6» 
Jo 

/■V 

F(v>,fc) = / (1 - fc 2 sin 2 ey 1/2 do. 

Jo 

Particularly, 

2 

Aq(<P,1) = -</>• 

7T 

At the critical point fc — > 1 or y — > z, we have 

(z 2 -y 2 )II(r,fc) = (z 2 -y 2 )^(fc) 

ttz 2 (z 2 -x 2 ) 1 / 2 ( 2 \ 

(z 2 -y 2 )Il(s,k) = (z 2 -y 2 )K(k) 

tt(1+lj)z 2 (z 2 -x 2 ) 1 ' 2 



(57) 



(58) 



+ 



and 



h 


= h 




- 






Ax 






1 




h 


1 




I 



2y/(l ~u)[{1-uj)x 2 -2uz 2 } 



1 Vr I, 



1- -<p s ,(59) 

7T ' 



(60) 



+ 



4^/(1 - w)[(l - w)a; 2 + 2loz 2 ] 



Let 



g = sinh 2 i? 1 < 1, |i?i|<|i?|. 



Then, at the critical temperature T C (H) 
4<Z 



/a 



z 2 = 4(1 + g ) 2 (l + 2g)/g, 



sm 



= sm 



1 + 9 
l + 2 9 



|2«-1| AT+? 



h 



8(1 + 9) 
1 + 4g 2 



(1 + 2q) 



2q + l V 1 + 2q 
'7 



2(l + g)|l -2q\(2q + 3) V 1 tt^ 



. (61) 



7 



and 



Ci 

c 2 
c 3 



where 

dRi 

~d~K 

dR 

~dK 



(q + l)(4q 2 



3) 



2q 



(dR dR 
coth2i?i— -i - coth2i?— 
oK oK 

-2(1 + q) I coth 2R t ^ - coth 2R^ 

( g + l)(2 g -l)(2g + 3)(2 g + l) 2 
2g(4g 2 + 1) 

/ dR dR 

x coth2i?i— - - coth2.R— 
V oK oK 



^ - (62) 



2( e -47^c s inh4iC c -7) 



27 



e -4 7 x c CO sh4^ c - 1 e -^ K " - 1 ' 
e -4 7 Kc ( sinh 4i f c _ 2 7 cosh 2 2if c ) + a sinh 4aif c 
2(sinh 2 2aif c + e - 4 ^ cosh 2 2if c ) 
-atanh2ai4T c (63) 



where K c = J/kT c (H). Combining these results, we ob- 
tain after some algebra the following expression for the 
energy at the critical temperature: 

1 dR 1 



U C (H)/J = -7- 2atanh2<xFsT c 

1 BR 
+ - ( 2-coth2tf)- 



coth 2R 1 

< 1 
9> 2" 



2sinh2i?i dK 
dR! 



dK 



coth 2i? 



dR\ 



In the limit of 7 — > for which R\ 
7r/4, we obtain 



dK) 1 
(64) 

00, </? r = </j s = 



U C (H)/J =-{2- V2)atanh2aK c 



7=0 



1 sinh AK C 



i sinh AaK r 



2-V2 sinh 2 2aiT c + cosh 2 2K C ^ 
which is the Fisher result jl|. Particularly, for H = 0, 
ASH gives the value -2\/l + V2- 



VI. MAGNETIC PROPERTIES 

For magnetic properties we take the derivative in 14511 
with respect to H and obtain 

M(if) = M C (H)+ const \T-T c (H)\ln\T-T c (H)\, 
T -► T c (iT) (66) 

where M C {H) is the magnetization at the critical tem- 
perature T C {H) given by 



M C (H) 



1 (yv±<^ s ) (l + e 4 ' K )tanh2a/f c 



< 1 
q> 2 



(67) 




1.5 2 2.5 3 3.5 
Temperature 1/K 

FIG. 6: Magnetization M versus temperature T at fixed mag 
netic field a for 7 = — 1. The broken line indicates M C (H). 




0.5 1 1.5 

Magnetic field a 



FIG. 7: Magnetization M versus magnetic field a at fixed 
temperature T for 7 = — 1. The broken line indicates M C (H). 

The magnetization M is plotted in Figs. 6 and 7 versus 
T and H, respectively, for 7 = — 1. To analyze M(H) = 
M( 7 , a, 1/K) at low temperatures, we note that 



M( 7) a,0) = 0, 
= 2, 



a < 1 — 7 
ev > 1 — 7 



(68) 



The value of M (7, a, 0) at a = 1 — 7 depends on how 
the zero temperature T = is approached. To see what 
happens we write 

a(K,0 = 1-7- ^ In 2(^ + 1) (69) 

where £ is a parameter controlling the approach to T — 
0. Particularly, the critical curve l|34fl is indicated by 
£ = 1. Along l(6"9"|) and K — > 00, one has i?i ~ 2JsT — > 00, 
tanh L — ► 1, and 



-4H 



1 + [2(V2 + !)]«, Jf-nx). 



8 



After some algebraic manipulation, we find 
m(£) = lim MUa(K,0,l/K 

K — >oc V 



„ "■■ - ;|) "1 + 1(1 -k)K(k) ,70) 

-I 7T 



with k = sinh 2e 2R. Here e = 1 for £ < 1 while e 
for £ > 1. We have 



-1 



m(£) = 

= 2-V2 
3 21 



£ = oo 
£ = 1 

4 ' 16*" (I)* 1 * 408 836 e = ° 
2 £=-oo.(71) 



-if i 




FIG. 8: Susceptibility as a function of temperature for 7 = 
-1. Solid line H = 0. Broken line H — J (a — 1/2). 



The susceptibility xCTja^^ -1 ^ which is the further 
derivative of (|flf))l with respect to H, diverges logarith- 
mically at T C (H) for H 7^ 0. The susceptibility x is 
plotted in Fig. 8 versus temperature T for 7 = — 1 and 
a = 0,1/2. The zero-field susceptibility x(7, 0, if" 1 ) is 
continuous. For example one has for 7 = — 1, 



X (-l, 0, K- 1 ) =2K-^{1- e 4R ) J 2 - I coth 2R 



+ coth 2i? 



2 cosh 2i?(3 - sinh 2 2R) + 6 



sinh 4R h 



(72) 



where e _2i? = (e AK + l)/2. The expression Ij72(l bears no 
direct relation with the critical energy (|53|l as found to 
exist in the case of 7 = 0. By direct differentiation 
of 16 7|) and making use of (|37|l . we obtain the critical 
susceptibility 



xhAK- x ) 



1 



l--((f r ± (f s ) 
7T 



< i. 



x (1 + e 47 ^ cosh" 2 2K C ), q > |73) 
Acknowledgment 

We thank R. Shrock for his interest and a critical read- 
ing of the manuscript. Work has been supported in part 
by NSF Grants PHY-0098801 (WTL) and DMR-9980440 
(FYW). 



[1] M. E. Fisher, Proc. Roy. Soc. A 254, 66-85 (1960); ibid. 

256 502-513 (1960). 
[2] H. Giacomini, J. Phys. A 21, L31 (1988). 
[3] K. Y. Lin, J. Phys. A 22, 3435 (1989). 
[4] K. Obredors, A. Labarta, A. Isaigue, J. Tejeda, J. Ro- 

drigues, and M. Pernet, Solid Stat. Commun. 65, 189-192 

(1988). 

[5] P. Azaria and H. Giacomini, J. Phys. A 21, L935 (1988). 
[6] C. Fan and F. Y. Wu, Phys. Rev. B 2, 723-733 (1970). 
[7] I. Syozi, Prog. Theor. Phys. 6, 306-308 (1951). 
[8] I. Syozi, in Phase Transitions and Critical Phenomena 

Vol. 1, Eds. C. Domb and M. S. Green (Academic Press, 

New York 1972). 



[9] R. M. F. Houtappel, Physica 16, 425-455 (1950). 

[10] Expressions of h,l2,l3 given just above Eq. (41) in |(| 
contain the following typos: The arguments g of the el- 
liptical integrals should read q. In the expression of I\ 
replace (z 2 — x 2 ) — » (z 2 — x 2 ) 1 ^ 2 , in the expression of I2 
multiply the whole expression by 1/z, and in the expres- 
sion of ^3 replace the first g by y. 

[11] M. Abramowitz and LA. Stegun, Handbook of mathemat- 
ical functions with formulas, graphs, and mathematical 
tables (National Bureau of Standards, Applied Mathe- 
matics Series 55, 1970), Eqs. (17.7.14) and (17.4.40). 



